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ABSTRACT 

Strong gravitational lensing magnifies the flux from distant galaxies, allowing us to detect 
emission lines that would otherwise fall below the detection threshold for medium-resolution 
spectroscopy. Here we present the detection of temperature-sensitive oxygen emission lines 
from three galaxies at 2 < z < 3.5, which enables us to directly determine the oxygen abun- 
dances and thereby double the number of galaxies at z > 2 for which this has been possible. 
The three galaxies have ~10% solar oxygen abundances in agreement with strong emission 
line diagnostics. Carbon and nitrogen ratios relative to oxygen are sub-solar as expected for 
young metal-poor galaxies. Two of the galaxies are Lya emitters with rest-frame equivalent 
widths of 20 A and 40 A, respectively, and their high magnification factors allow us for the 
first time to gain insight into the physical characteristics of high-redshift Lya emitters. Us- 
ing constraints from the physical properties of the galaxies, we accurately reproduce their 
line profiles with radiative transfer models. The models show a relatively small outflow in 
agreement with the observed small velocity offsets between nebular emission and interstellar 
absorption lines. 
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1 INTRODUCTION 

Strong gravitational lensing of distant galaxies behind massive 
galaxy clusters allows us to study intrinsically fainter galaxies 
than possible for regular field galaxies. With magnification factors 
reaching up to 50, it becomes possible to study the detailed phys- 
ical properties of galaxies that are fainter than the characteristic 
luminosity, L* , at a given redshift. Since the first de tailed analy- 
ses o f a lensed Lyman break galaxy, MS 15 12-cB58 ( Pettini et al.l 
l200d : iTeplitz et alj |2000| : IPettini et alj l2002h . a steady increase 
in the number of spectroscopic observations of bright galaxies 
at z ~ 2 either detected se rendipitously or in l arge dedicated 
surveys have appeared (e.g. iFosburv et al.1 120031 : ICabanac et alj 



* Based on data from the X-shooter GTO observations collected at the 
European Southern Observatory VLT/Kuyuen telescope, Paranal, Chile, 
collected under programme IDs: 084.B-0351(D), 086.A-0674(A), 086.A- 
0674(B), 087A-0432(A) and 087A-0432(B). 
f lise@dark-cosmology.dk 



20051: ISwinbank et al.l [20071: lAllam et alj|2007l: [stark et alj|2008l: 
Lin et alj|2009l: iDiehl et al.ll2009l : lHainline et aI1l2009l : iBian et al.1 
20ld:lRigbvetal.ll201l|) . To date, most lensed galaxies that have 



received more attention still belong to the bright end of the lu- 
minosity function, although recent investigations h ave started to 
probe the lower-mass en d of the galaxy distribution dRichard et al.1 
1 20 111 : IWuvts et alj|2012h . Some galaxies show particularly strong 
UV emission lines, suggest ing the presence of a large population 
of ve r y massive, hot sta rs dFosburv et alJl20ol&IVillar-Martm et al.1 
2004; iRaiter et al. 2010), giving rise to nebular emission lines that 
are very rarely detected in nearby galaxies. Other strongly star- 
forming lensed galaxies have UV spect ra that do not differ sub- 
stantially from the cB58 galaxy spectru m dOuider et al.l200sj|2010l : 



stantially trom the cB3S galaxy spectrur 
iDessauges-Zavadskv et al.ll2010Ll201 lh 



With the boost of the fluxes from gravitational lensing, it also 
becomes feasible to detect emission lines such as [Olll] A4363 
or [N III] A5755 which, when related to other transitions from the 
same ionisation stage, are sensitive to the gas temperature. Com- 
bining this with the knowledge of the density of the gas, it is 
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possible to derive the abu ndances of various chemical elements 
directly dOsterbrockl [l989). This approach is widely used to de- 
termine abundances in low-redshift galaxies. The flux from the 
[O III] A4363 line increases with decreasing metallicity; however, 
low-metallicity galaxies at high redshift are also very faint, so the 
line is difficult to detect, in particular since the line is redshifted 
to near-IR part of the spectrum. The high-redshift L yman break 
galaxies (LBGs) with <25.5 mag dSteidel et al.tf l996)) are known 
to be massive galaxies, and because of the relati on between mas s 
and metallicity, LBGs are typically metal rich ( lErbetal.ll2003) . 
This implies that their temperature-sensitive lines have fluxes be- 
low the detection threshold even on the largest telescopes. To date, 
only one dete ction of [Qui] A4363 from a galaxy at z > 1 has 
been reported dYuan & Kewlevl2009b . while ot her attempts to mea - 
sure the line have led to upper limits only dRigbv et alj 1201 lh . 
Analysing other rest-frame UV lines (O III] AA1661,1666), which 
are also sensitive to the gas temperature may alleviate this difficulty 
of detecting [Qui] A4363 and determini ng abundances directly 
dVillar-Martm et alj|2004l ; lErb et alj[201oh . However, this doublet 
is also rarely observed in high-redshift galaxies, and it is more sen- 
sitive to uncertainties in the correction of extinction in the UV range 
of the spectrum. 

To avoid the need for very long observation times to detect the 
temperature sensitive lines, alternative methods have been devel- 
oped. Using only strong rest-frame optical emission lines, various 
diagnostics have been cali brated with respect to the direct abun- 
dance determina tions (e.g., iPagel et alJl979l : |Pettini & Pagell2004l ; 
iKewlev & Ellisonll2008l) . Since the temperature-sensitive lines are 
generally too weak to be measured even in bright emission-line 
galaxies at high redshifts, we are forced to use the relations estab- 
lished at low redshifts to derive abundances. The validity of using 
the exact same strong-line diagnostics at high redshifts has to be 
verified. 

To find potential UV bright high-redshift galaxies, we have 
studied a sample of gravitationally lensed galaxies at 1 < z < 6, 
based on selecting high surface brightness regions in extended arcs. 
In a companion paper dChristensen et al.ll2012L hereafter Paper I), 
we present the sample and the continuum emission properties of 
the galaxies. We determine stellar masses from spectral fitting, star 
formation rates from emission lines, and oxygen abundances. We 
also demonstrate that the low-mass galaxies exhibit a large scat- 
ter compared to the fundamental relation for st ar-forming galaxies 
dLara-Lopez et alj|2010l : iMannucci et alJlzblOh and provide an al- 
ternative calibration valid for low-mass galaxies. 

Out of the 12 galaxies in our original sample, three galaxies 
at 1.8 < z < 3.5 show temperature-sensitive lines. Two of these 
are also strong Lya emitters, while the third galaxy shows Lya in 
absorption. In this paper we present a more detailed analysis of the 
nebular emission lines from these three galaxies. We present the 
data in Section [2] and derive the direct abundances and determine 
escape fractions in Section|3] In Section|4]we model the two Lya 
emission lines using constraints from the physical properties of the 
galaxies. 



2 OBSERVATIONAL DATA 

The data for the lensed galaxies were obtained with X-shooter on 
the VLT during separate obser ving runs in 2010-201 1 . One of the 
scien ce drivers for X-shooter dD'Odorico et aljfeood IVernet et al.l 
l201ll) is to determine redshifts of faint single targets and anal- 
yse emission line galaxies at z > 1. In previous papers, 



we demonstrated t he ability of X-shooter for the observations 
of lensed galaxies 1 Dessauges-Zavadskv et al.1 l2010l : IPettini et al.l 
l20ld : IChristensen et alj|20ld) . In Paper I, we describe the target 
selection, observations, data reduction, and analysis of the spectra. 
In this paper, we focus on the analysis of the final emission line 
fluxes from three galaxies after correcting for many effects: Galac- 
tic extinction, intrinsic extinction, slit losses, and lens magnifica- 
tion. One of the lensed galaxies is found in the well-studied Abell 
1689 cluster, while the remaining two are found in tw o southern 
equivalents to the MACS clusters dEbeling et alfeoOll) with decli- 
nations < - 40 deg (H. Ebeling, unpublished). This section presents 
the three galaxies which we can investigate in further detail. 

2.1 Abell 1689 arc ID 31.1 



The arc ID 31.1 is adopted from iLimousin et all d2007h . For this 
target, a wealth of emission lines are observed as listed in Ta- 
ble [TJ From the strongest emission lines, which are not affected 
by telluric absorption, and also excluding Lya we find the redshift 
z = 1.8339 ± 0.0003. The last column in Table Q] lists the line 
fluxes after correcting for intrinsic extinction as described in Sec- 
tionl3~T1 

The Lya line profile is double-peaked with a much brighter 
red than blue component, which is the characteristic signature of 
resonance scattering of photons in an outflowing medium (see Fig. 
|2j- Its rest-frame equivalent width is 40A. In Sect. [4] we model the 
line profile taking into account the constraints from physical pa- 
rameters determined from other emission lines. 

We detect the temperature-sensitive [O III] A4363 line. Its 
strength increases relative to [O III] A5007 with decreasing metal- 
licity, but at high redshifts, the galaxies that have the best signal- 
to-noise ratio spectra are also the brightest and most massive, 
and hence most metal-rich ones at those redshifts. To date, [O III] 
A4363 has only been observed in a nother lensed galaxy at z = 1.7, 
also behind the Abell 1689 cluster dYuan & Kewlevll2009l) . 

Many other rest-frame UV emission lines are identified in the 
Abell 1689 ID 31.1 spectrum. In particular we draw attention to the 
narrow emission lines from the Civ AA1548,1550 doublet shown 
in Fig. [TJ Civ AA1548,1550 are almost always seen in absorp- 
ti on in high-redshift galaxies (another ex ception is the Lynx arc 
in lHolden et"al]|200ll ; ¥osburv et alj|2003l) . or in P-Cygni profiles 
with a minor cont ribution from the em ission component compared 
to the absorption dPettini et al]|2000l) . This indicates a highly un- 
obscured line of sight to the galaxy. Searching for the UV absorp- 
tion lines in the low signal-to-noise ratio spectrum, we only de- 
tect weakly the Sill A1260,1264, and 1526 lines at z = 1.8348 ± 
0.0012, while we do not detect the typically strong high-ionisation 
absorption doublet Si IV A1393,1402. Due to the weak absorption 
lines, we argue that the interstellar medium (ISM) only has a par- 
tial coverage, like observed for a couple of other lens ed galax- 
ies dOuider et al. 2009; Dessauges-Zavads kv et al. 2010). Among 
the uncommon emission lines, we also detect the temperature- 
sensitive Olll] AA1661,1666 lines, which among other high-z 
galaxies are also observ ed in the bright z = 2.3 galaxy Q2343- 
BX418 terb et alJl201Ch and in the Lynx arc dFosburv et alj|20()3 : 
Ivillar-Martm et al.ll2004l) . 

2.2 SMACS J0304.3-4402 ID 1. 1 

The X-shooter spectrum covers one of the multiple lensed images, 
and within the slit three distinct regions are visible in the HST im- 
age (see Paper I). The emission lines listed in Table|2]present a sum 
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Table 2. Emission lines from the SMACS J0304 arc 



4370 



4380 



4390 
Wavelength (A) 



4400 



4410 



Figure 1. Section of the observed UVB spectrum of the ID 31.1 arc show- 
ing its highly uncommon narrow Civ AA1548.1550 emission lines. This 
suggest an ISM with a low covering factor along the line of sight. 



Table 1. Emission lines from the A 1689 31.1 arc 



line 




Aobs 


f\ 

J obs 


f b 

J cor 


Lyce 


1215.67 


3445.77 


72.7±0.7 


4309±42 


CIV 


1548.20 


4386.81 


4.3±0.4 


137±13 


Civ 


1550.77 


4394.22 


4.7±0.3 


149±10 


Oiii] 


1660.81 


4705.82 


3.2±0.4 


90±11 


Om] 


1666.15 


4720.63 


7.4±0.3 


206±8 


[Cm] 


1906.68 


5402.29 


9.1 ±0.4 


202±9 


Cm] 


1908.73 


5408.19 


4.7±0.4 


104±9 


[On] 


3727.09 


10563.06 


9.2±0.8 


68±6 


[On] 


3729.88 


10569.40 


11.7±o.9 


86±7 


[Ne in] 


3869.84 


10965.28 


6.4±2.5 


44±17 


H7 


4341.69 


12301.90 


17.0±2.1 


97±12 


[Om] 


4364.44 


12366.36 


6.2±1.8 


35±10 


H/3 


4862.70 


13777.99 d 


42.7±4.9 


205±24 


[Om] 


4960.29 


14057.55 d 


63.0±10.6 


293±49 


[Om] 


5008.24 


14194.09 c 


211.5±5.7 


968±26 



Notes. 

a Rest-frame vacuum wavelengths. 

b Emission line flux in units of 10~ 18 erg cm -2 s~ x . / b s and / cor cor- 
respond to the observed flux and that corrected for intrinsic absorption, re- 
spectively. 
c Close to sky line. 

d Affected by strong telluric absorption lines. 

of the total flux from the components within the slit. The redshift 
for this galaxy is measured to be 2 = 1.9634 ± 0.0002. Again, we 
detect the O III] AA1661,1666 doublet. Apart from velocity offsets 
of ~100 km s _1 between emission lines from the two components, 
there are no significant differences in the relative line fluxes as ver- 
ified by comparing extracted spectra from the two distinct compo- 
nent. The UV continuum has the same flux, but redwards of 2600 A 
in the rest-frame, the continuum emission differs by a factor of two, 
with the smallest flux coming from the most compact component. 
The emission line fluxes listed in Table [2] are the total flux in the 
X-shooter slit summed for the two components. Even though the 
Balmer lines are very bright for this source, Lya emission is ab- 
sent, and its (UV ) spectrum appears similar to the cB58 spectrum 
jPettini et al]|2000h . 

In addition to the strong UV absorption lines in the ISM, a 
strong intervening Mg II at z = 1.5492 is seen in the VIS spectrum 



line 


A r a 


A ODS 


f b 

■> obs 


/cor 


Olll] 


1660.81 


4921.93 


7.1±1.6 


56±13 


Om] 


1666.15 


4937.35 


11.4±1.9 


90±15 


[On] 


2470.22 


7322.39 


8.4±2.1 


45±11 


[On] 


3727.09 


11044.78 


437.8±3.4 


1514±12 


[On] 


3729.88 


11053.09 


587.2±3.9 


2029±14 


H9 


3836.49 


11370.55 c 


29.1 ±4.2 


98±14 


[Ne in] 


3869.84 


11468.93 c 


149.9±3.8 


499±13 


He 1 + H8 


3890.17 


11527.64° 


95.3±3.2 


315±11 


[Nem] 


3968.53 


11761.36 


64.3±3.4 


209±11 


H7 


3971.20 


11767.48 


51.8±1.7 


168±6 


H<5 


4102.92 


12158.43 


128.5±2.7 


403±9 


H 7 


4341.69 


12866.19 


205.9±1.3 


610±4 


H/3 


4862.70 


14410.72 c 


500.0±0.9 


1324±2 


[Om] 


4960.29 


14700.41 c 


669.5±1.9 


1739±5 


[Om] 


5008.24 


14841. 33 c 


2357±3.5 


6066±9 


He I 


5877.28 


17417.44 


47.0±5.2 


104±12 


Ha 


6564.63 


19453.85 c 


1802.1 ±0.6 


3645±1 


[Nil] 


6585.42 


19513.96 c 


58.1±2.7 


117±5 


[SII] 


6718.29 


19910.44 


156.8±1.1 


311±2 


[Sll] 


6732.67 


19951.76 


102.7±1.7 


203±3 



Notes. 

a Rest-frame vacuum wavelengths. 

b Emission line flux in units of 10~ 18 erg cm -2 s _1 . / ODS and / cor 
correspond to the observed flux and that corrected for intrinsic absorption, 
respectively. 

c Affected by telluric absorption lines. 



where its lines are blended with the strong Fe II lines of the lensed 
galaxy's ISM. 



2.3 SMACSJ2031.8-4036ID1.1 

Emission lines from this arc are listed in Table [3] The detected 
lines, excluding Lya, give a flux-weighted average redshift of 
z = 3.5073±0.0002. At the low S/N per pixel in the UVB and VIS 
arms, we can only identify several Lyman forest absorption lines, 
and Sin A1206 and Siv AA1394,1403 at z = 3.5061 ± 0.0013. A 
detailed analysis of the absorption lines requires a spectrum with 
a higher S/N measured in the continuum, so we refrain from any 
detailed analysis of the present data. Again, the Lya emission line 
profile is double-peaked, with a rest-frame EW of 20 A (see Fig.[2]l, 
and we will return in Sect. |4]with a more detailed treatment of this 
line. 

At a low S/N we detect the semi-forbidden Niv] A 1486 line 
which is rarely seen in high-z galaxies, although a few detections 
of this line in emission ar e reported for lensed galaxie s, for example 
cB58 jPettini et ai1l2000h , the Lynx arc at z = 3.3 jFosburv et all 
120031) . and a galaxy at z = 5.563 dRaiter et al.ll2010h . When this 
line is detected, it can be interpreted as being caused by a large 
population of young and very hot, massive stars, but also an ob- 
scured AGN could be responsible for the ionisation (Vanze lla et al.l 
l2010bh . Since neither Nv AA1238,1242 nor Hell A1640 are de- 
tected, we exclude the po ssibility that the ion isation is caused by 
an AGN (see Table 1 in iBinette etakl 120031) . Hidden AGN like 
in radio galaxies typical ly display strong Hell A1640 lines (e.g. 
Ivillar-Martm et al.lll999l) . 
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Table 3. Emission lines from the SMACS J2031 arc 



line 


A,. a 


-^obs 


fb 

■> obs 


/cor 


Lya 


1215.67 


5480.80 


173.7±0.5 


242±1 


NIV] 


1486.50 


6700.85 


7.6±0.9 


10.1±1.2 


O III] 


1 660 8 1 


7485.15 


2.9±0.6 


3.8±0.8 


Oin] 


1666.15 


7509.46 


8.8±0.7 


11.5±0.8 


[Cm] 


1906.68 


8594.04 


12.0±0.5 


15.4±0.6 


Cm] 


1908.73 


8603.54 


8.4±0.9 


10.8±1.2 


[On] 


3727.09 


16800.18 


14.4±1.6 


16.9±1.9 


[On] 


3729.88 


16811.01 


15.5±1.3 


18.2±1.5 


[Ne III] 


3869.84 


17442.96 


14.8±2.1 


17.3±2.5 


Hl(H7) 


3971.20 


17898.76 


6.0±1.0 


7.0±1.2 


H/3 


4862.70 


21917.47 


37.7±1.3 


42.8±1.5 


[Om] 


4960.29 


22357.68 


61.9±0.9 


70.1±1.0 


[Oin] 


5008.24 


22573.68 


205.2±0.5 


232±0.6 



Notes. 

a Rest-frame vacuum wavelengths. 

b Emission line flux in units of 10 — 18 ergcm~ 2 s _1 . / t, s and / cor 
correspond to the observed flux and that corrected for intrinsic absorption, 
respectively. 



3 PHYSICAL CONDITIONS FROM EMISSION LINES 

All emission line fluxes listed in Tables 1-3 are correc t ed for Galac- 
tic extinction using the dus t maps in [ Schlegel et al. ( 1998) along 
with an extinction curve in iFitzpatricld H999h with Rv=3.l. All 
the physical properties derived in this section are summarised in 
Table[4] When line ratios are considered, we do not apply any cor- 
rection for slit losses and magnification factors, since we assume 
that the emission lines come from the exact same region. 



3.1 Gas phase reddening 

Any ratio of Balmer lines can be used to derive the reddening in the 
gas phase of the galaxies. In the absence of reddening, the expected 
emission-line ratio is a function of the gas temperature and den- 
sity. We adopt the values tabulated for t emperatures T = 10 000- 
20 000 K and densities n ~ 100 cm" 3 in lBrocklehurs3 dl97 lh . The 
choices of these values are based on the estimated gas temperature 
and density in Sect |3.3| 

To derive the gas phase reddening, E(B — V) gas , we use all of 
the Balmer line pairs available and calculate the weighted average, 
such that the reddening is determined by the Balmer line ratio with 
the smallest uncer tainty. We assume an extinction curve derived for 
starburst galaxies (Cal zettiet alj2 000). and the resulting reddening 
values are listed in Table [4] The reddening determined for M2031 
is uncertain, and a value of zero is consistent with that determined 
from fitting the entire spectrum with stellar population models as 
demonstrated in Paper I. 



3.2 Starburst/ AGN ionisation 

The ionising radiation from either AGN or recent formation of mas- 
sive stars gives rise to strong emission lines in galaxies. Whether 
the narrow emission lines are caused by an AGN rather than a 
strong starburst is usually investigated through emission-line ra- 
tios, which are sensitive to the hard ionisation from AGN. When 
only the strongest rest-frame optical emission lines can be detected, 
as is typically the case in high-redshift galaxies, the line ratios 



[Om]/H/3 versus [N Il]/Hq can be used to distinguish the domi- 
nant ionising source JKewlev & Dopita 2002). Among the galaxies 
anaysed in this paper, only M0304 have all the relevant lines de- 
tected, and its line ratios are characteristic of star-forming galaxies. 

In M2031, the C IV AA1548,1550 lines are not detected 
to a level of 3 x 10~ 18 erg cm -2 s -1 , so the ratio Civ 
AA1548,1550/[Ciii]Ciii] AA1907,1909<0.15 indicates a softer 
ionising spectrum, w hile an AGN would cause a ratio of 2 as ar- 
gued for the Lynx arc dBinette et al.l2003T) . The line ratio for A3 1 . 1 
is Civ/Cm] = 0.65±0.71 indicating a higher ionisation parame- 
ter. In AGNs, the lines N V AA1238,1242 and He II A1640 are ex- 
pected to be strong, and since these lines are not detected in any of 
our spectra it suggests that the main contributor to the ionisation is 
massive stars. 

3.3 Gas temperatures and densities 

Although [O III] A4363, which is conventionally used for direct 
temperature measurements in low-redshift galaxies, is not de- 
tected in the M2031 source, the electron temperature can instead 
be derived directly from the Om] AA1 661,1666/[Q hi] A5007 ra - 
tio as demonstrated for the Lynx arc {Villar-Martm et al. 2004). 
We use the nebular package in IRAF JShaw & Dufourll 19951) 
to determine the gas phase temperatures and densities. From the 
oxygen line-flux ratio, we find T e — 16100 ± 300 K. The 
[C III],C III] AA1907,1909 doublet can be used to derive densi- 
ties. For M2031 we determine n e = 2900 ± 5000 cm" 3 given 
that T e = 16100 K. For the M0304 source, we also use the 
O III] AA1661,1666 / [O III] A5007 ratio to derive the temperature 
and density directly, again after correcting the measured fluxes for 
reddening. The results are listed in Table [4] 

Two different diagnostics can be used to derive the elec- 
tron temperature for the Abell 1689 31.1 source. Using the 
Om] AA1661J666 lines gives T e = 26900 ± 1000 K, while the 
commonly used [Om] 4363 gives T e = 2 1600 ± 300 K as 
determined iteratively using the equations in lAllerl dl984h . Since 
the ratio from [O III] A4363 over [O III] AA4959,5007 is less af- 
fected by reddening, we chose to use the latter temperature mea- 
surement for further analysis of the oxygen abundance. The ratio 
of the [C III],C III] AA1907,1909 (1.9 ± 0.2) doublet is unphysical 
at the 1.2cr level, as the maximum allowed fraction is 1.65, how- 
ever, the high value is consistent with a low electron density. 

Since the [O II] AA3727,3730 doublet is detected in all three 
galaxies, we can also use this to determine the electron density as 
listed in Table [4] The accurately measured electron densities from 
the medium resolution data of the three galaxies is similar to the 
densities in other lensed z ~ 2 galaxies whe re a range between a 
few 100 and a few 1000 cm" 3 are common faaimineetail 2009; 
IWuvts et al.ll2012h . 

3.4 Abundances 

In this section we derive the integrated abundances of oxygen and 
carbon. These values are not dependent on magnification, because 
these parameters are derived from line-flux ratios rather than abso- 
lute values. 

3.4.1 Direct (T e ) abundance measurements 

Since we have measured the temperature directly for the sources, 
we can determine the oxygen abundance with direct temperature 
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methods. We use the equ ations in llzotov et alj l l2006h with the 
atomic data referenced in lStasinskH J2005h to derive 2+ /H + and 
+ /H + . The sum of these two contributions gives the total oxy- 
gen abundance, when neglecting contributions from higher ionised 
stages, as is conventionally done. We take into account that the two 
different ions may reside in different zones of the H II region, and 
therefore have different electron temperatures. Since we do not de- 
tect the [Oil] AA7320,7331 lines in any of the spectra, we assume 
a relation betwee n the tempera t ures T (Q II) and T(0 III) described 
by equation 14 in llzotov et al.l d2006h . The resulting oxygen abun- 
dances are listed in Table[4] 

Using the temperature and densities determined in the section 
above, we calculate the total oxygen abundance for M2031. The 
derived uncertainty includes both the uncertainty from the emis- 
sion line fluxes, temperature, and density. The contribution from 
singly ionised oxygen is 12+log(0 + /H + ) = 6.94±0.06 and double 
ionised oxygen is 12+log(0 2+ /H+) = 7.69 ± 0.02. The sum gives 
the total oxygen abundance 12+log(0/H) = 7.76 ± 0.03, assuming 
that other oxygen ionisations states do not contribute. 

For M0304 the contribution from single and double- 
ionised oxygen is 12+log(0+/H + ) = 7.61 ± 0.11 and 
12+log(0 2+ /H+) = 7.89 ± 0.08, respectively. The sum gives 
the total oxygen abundance 12+log(0/H) = 8.07 ± 0.09, again 
assuming that other oxygen ionisation states do not contribute. 

For Abell 1689 ID 31.1 the contribution from single and 
double-ionised oxygen is 12+log(0+/H+) = 6.42 ± 0.33 and 
12+log(0 2+ /H + ) = 7.66 ± 0.11, respectively, and the total oxy- 
gen abundance is 12+log(0/H) = 7.69 ± 0.13 with the TXOm) 
temperature determined from [O III] A4363 in Sect. |3.3l 

3.4.2 Strong-line diagnostics 

In high-redshift galaxies temperature-sensitive lines are typ- 
ically not detected. Instead one has to rely on the rela- 
tions between strong emission line ratios and the directly 
derived oxygen abundances. The ratio defined as R23 = 
([OII1AA3727, 37 30 + [OIII]AA4959, 5007) /H/3 introduced by 
|Pageletal.l ( fl97^) is widely us ed, and its calibrat i on relies either on 
photo-ionisation models (e.g. , lMcGaugrifl99lU Kobul nickv et al.1 
1 19991 : iKewlev & Dopit3 [2002). or comparison with objects where 
direct oxy gen abundan ces are derived from te mperat ure- sensitive 
lines (e.gjAlloin et al.lll979l;|Pilvugiir& Thuan 2005; Naga o et all 



2006). The R23 relation has two possible solutions with a 
high-metallicity and a low-metallicity branch. A degenerate so- 
lution can be avoided if other emission lines, such as [Nil] 
A6586 or [Nem] A3869, are detected dDenicolo et alj |2002| ; 
IKewlev & Dopitall2002l ; IPettini & Pageill2004 iLiang et alj|2006h . 
Also the [O ill] A5007/[0 11] A3727 line ratio can be used to dis- 
tingui sh between th e upper and lower branch of the R23 calibra- 
tions ( Nag ao et al J I200Q) . Any strong emission line diagnostics 
have large intrinsic scatters of ~0.2 dex, while v arious calibrations 
may differ by up to 0.7 dex relative to each other ( Kew lev & Ellisor] 
2008), thereby making it difficult to compare objects from the lit- 
erature. 

The determinations of oxygen abundances using strong-line 
diagnostics were presented in Paper I, where we used [O II], [O III], 
or [Nelll] emission line ratios to distinguish between the up- 
per and lower branch of the R23 calib ration using the met hod in 
IPilvugin & Thuanl d2005h, the 03N2 in lPettini & Pagell d2004l) and 
Ne3Q2 in lNagao et al.l d2006l) . For direct comparison, the derived 
oxygen abundances are reproduced in Table [4] In addition to the 
uncertainty of the abundance derived from the emission line ra- 



tios, the calibrations themselves have an intrinsic scatter of ~0.2 
dex, which must be taken into account when comparing the strong- 
line diagnostics with the direct measurements. In any case, it is 
encouraging that the stron g-line calibration, in particular the one in 
IPilvugin & ThuarJ 12005]) agrees very well with the direct oxygen 
abundances. 



3.4.3 C/O abundance 

The abundance ratio C/O increases with increasing oxygen abun- 
dance above 10% solar. Th is trend is detected in halo stars, H II re- 
gions dGarnettetalll 19951). an d in th e integrated spectra of Lyman 
break galaxies dShaplev et al.ll2003h . A sub-solar C/O ratio arises 
when the production of O and C is dominated by SN Type II ex- 
plosions, while C can also be produced by intermediate mass stars. 
Another explanation for the variation of the C/O ratio as a function 
of the oxygen abundance involves a metallicity depend ence of the 
mass-loss rate for high-mass stars 1 Akerman et al.|[2004h . 

Following the method in lErb et alj f2oToh we use the ratio of 
Oin] AA166U666 and [Ciii],Cm] AA1907,1909 emission lines 
to determine the C/O abundance ratios. Since the ionisation poten- 
tials for oxygen is higher than for carbon, C 3+ could be present in 
a medium with a high ionisation parameter. The ionisation parame- 
ter U is defined as the ratio between the density of ionising photons 
and the density of hydrogen atoms. 

In the spectrum of M0304, the [C III] ,C III] AA1907,1909 dou- 
blet is not detected, even though it could have been, if the lines 
were as bright or brigh ter than O Hi] AA166 1,1666. This is con- 
sistent with the note bv lShaplev et alj d2003l) that the [C III],C III] 
lines appear to be stronger in LBGs with Lya in emission. 

Once we know the metallicity of a galaxy, the ratio of the 
emission lines [Om] AA4959,5007/[O II] AA3727,3730 can be 
used to determine the ionisation para meter using the photoioni - 
sation models and iterative equations in IKewlev & Dopital d2002l) . 
The three galaxies all have logf/ between -3 and -2 as listed 
in Table [4] Similar high ionisation paramet ers have been mea- 
sured for other grav i tationally lensed g alaxies dHainline et alj2009l 
iRichard et alj|201 luRigbv et aU201 II) . as well as changes in other 
emission lines ratios in LBGs relative to l ocal galaxies have be en 
attributed to higher ionisation parameters (Erb et al. 2006, 2010). 

For an ionisation parameter of log U « —2, the dominant 
ionisation level for both C an d O are d ouble ionised, and one can 
approximate C/0~C 2+ /0 2+ dErb et alj2010t) . while at higher ion- 
isation parameters, an ionisation correction factor (ICF) is neces- 
sary to include. Given the values of the ionisation parameters in 
the three lensed galaxies, no ICFs are needed to derive the C/O ra- 
tio. For M2031 we calculate log(C 2+ /0 2+ ) = -0.80 ± 0.09, while 
for Abell 1869 ID 31.1 we derive log(C 2+ /0 2+ ) = -1.03 ± 0.08 
based on the same emission line ratio diagnostics. The results are 
listed in Tab le [4] Compared to the solar value log(C/O)0 = -0.26 
( Asplund et alj 20091) , the observed ratios are sub-solar. A decreas- 
ing C/O ratio wi th decreasing metal l icity is a know n property of 
un-lensed LBGs dShaplev et al.ll2003l : lErb et alfeoiol) . and the low 
C/O ratios of the lensed galaxies measured here are consistent with 
their ~0. 1 solar oxygen abundance. 

The fact that we observe a low C/O ratio when the galaxy 
ages are very young (~4 Myr according to SED models in Pa- 
per I), suggests that intermediate-mass stars h ave not contribut ed 
to the ISM in the galaxies. As concluded by lErb et alj d2010h a 
C/O value consistent with that observed in local metal-poor Hll 
regions suggests that the carbon production in high-mass stars de- 
pends on their metallicities as predicted by either nucleosynthe- 
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sis mo dels (Weaver & Wooslevll 19931) or by mass loss from stellar dening, similar to what we find for the two Lya emitters in this 
winds lMa ederill992 ). study. 



3.4.4 N/O abundance 



The nitrogen to oxygen ratio is an interesting parameter, because 
nitrogen provide s a measure of the c hemical enrichment history 
of a galaxy (see iPetitiean et al.l 12008). Nitrogen can be produces 
in two ways: a primary production from newly synthesised carbon 
via the CNO cycle, and a secondary production also from a CNO 
cycle, but with enhancement from a previous generation of stars 
that create additional carbon and oxygen. 

To determine the N/O abundance ratio for M0304, which is 
the only source in the sample where [N il] A6586 i s dete cted, 
we use the calibration of log(N+/0+) in iPagel et ail dl992l) . To 
account for the different ionisation stages arising in different 
zones in the H II region, we convert the tempe rature measured 
for T(OlIl) to T(Ol l) as in llzotov et all (2006). As argued by 
iThurston et alj dl996t) . ionisation models of low-metallicity gas 
show that log(N/0) ~ log(N + /0 + ), and that an ionisation cor- 
rection factor is not necessary. For the M0304 source we find 
log(N/0) = —1.64 ± 0.05, which is below the plateau seen 
in low-metallicity Hll r egions in dwarf and irregular galaxies 
dvan Zee & Hav nes 2006). The observed (N/O) ratio is more sim- 
ilar to the ratios observed in the more metal-rich (~0.1 so- 
lar) damped Lyman- a s ystems (Petitjean et al. 2008; Petti ni etHEI 
l2008l : lBattisti et alj2012 | ). and considerably below the solar level of 
log(N/Q)^ = -0.86 dAsplund et alj |2009). This sub-solar (N/O) 
ratio suggests that the galaxy, like in high-redshift damped Lyman- 
a absorbers, has only experienced the primary production of nitro- 
gen from the CNO cycle, and cannot have had a large population 
of stars forming earlier because these would have increased the ni- 
trogen level. This is consistent with a young age of ~60 Myr based 
on stellar population model fits in Paper I. 



3.5 Lyman-a escape fraction 

Two sources, M2031 and A 1689 ID 31.1 have double-peaked Lya 
emission lines, whic h is not uncomm on as it is seen in 30% of 
LBGs at z = 2-3 dKuIas et al.1120121) . and up to 50% of Lyman- 
q emitters dYamada et alj|2012r) . Only few le nsed galaxies studied 
in detail to date ar e also strong Lya em itters dFosburv et al.ll2003l ; 
lOuider et alj|2009l : iBavliss et aljfeoich . so the present data allow 
for a more in depth analysis including constraints from the rest- 
frame optical emission lines. The escape fraction can be determined 
from the observed Balmer emission lines relative to the observed 
Lya line flux. Under the assumption of a case B recombination 
scenario, and zero extinction, we find that the escape fraction of 
Lya photons from M2031 is 20%, and 4% or 85% for the A1689 
31.1 arc before and after correction for an intrinsic reddening of 
E{B-V) =0.2,7. 

At comparable redshifts, the comparison of galaxies UV lu- 
minosity to their Lya emission lines have revealed a median Lya 
escape fr action of 29%, wh ere any value between and 100 % is 
possible dBlanc et alj|201lh . Similarly to our calculations, compar- 
isons of the flux f rom Ha and Lya lines for both z ~ 0.3 gala xies 
dAtek et alj 2009) and a few z ~ 2 galaxies dHaves et al.ll2010h in- 
dicate a large range of Lya escape fractions from to 100%. There 
is a tendency that large r escape fractions are found for gal axies with 
smaller reddenings dHaves e t al. 2010; Bla nc et aU201 lh . but only 
when the Lya emission lines are not corrected for the intrinsic red- 



3.6 Lyman limit escape fraction 

The M203 1 source has a sufficiently high redshift that we can de- 
rive an upper limit on the escape fraction bluewards of the Lyman 
limit. In the region 3800 — 4100 A, corresponding to a rest-frame 
wavelength of ~ 900 A, the average flux measured in the extracted 
ID spectrum after binning by a factor of 50 in the dispersion is 
(1 ± 5) x 10~ 20 erg cm -2 s _1 A~ x , while the observed flux at 
the rest-frame 1500 A is (2.3 ± 0.3) x 10~ 18 erg cm" 2 s" 1 A -1 . 
The relative escape fraction at these wavelengths needs to be cor- 
rected for the absorption in the interg alactic medium (IGM) (e.g., 
Ilnoue et al]2006l ; lBoutsia et alj201 lh : 

, (fl50o/-f90o)int , IGMi 

/csc,ll = 77 77 r exp(-T 900 ), 

( J 1 500 //900 Jobs 

where T — exp(— Tg§) M ) is the transmission in the IGM at 
the rest-frame 900 A. The intrinsic fraction of the flux densities 
(Ji50o/-f90o)int is model dependent, and we determine the ratio 
from the best fit spectral template in Paper I to be 1.44. Since we 
do not find evidence for significant intrinsic reddening, this value is 
not corrected for extinction. From the observed ratio and an optica l 
depth of the IGM at z = 3.5 due to the Lyman series dMadauh995T) 
giving a transmission of 0.2 at A = 900(1 + 3.5) A, we derive an 
upper limit for the escape fraction / OS c,ll < 0.031 ± 0.025, or 
<3.0±2.5%. However, the IGM transmission varies significantly 
along different sight lines due to the random occurrence of inter- 
vening Lyman limit systems, and at z ~ 3.5 the IGM transmission 
just bluewards of the Lyman limit is T = 0.3 ± 0.2 based on nu- 
merical simulations dlnoue & Iwatai r2008). Accordingly, our limit 
to the escape fraction is less constraine d; / eS c,LL <11%. 

In a sample of 14 LBGs at z ~ 3 lShaplev et all d2006t) found 
a significant escape fraction in two galaxies implying an average 
escape fraction of 14%. However, because spatial offsets between 
the galaxy emission in the rest-frame Lyman continuum and at 
rest-frame optical wavelenghts have been measured, the large es- 
cape fraction c ould be contaminated b y lower redshift interlopers. 
Consequently, IVanzella et alj d2010al) find that LBGs at z ~ 4 
have upper limit of / csc ,ll < 5-20%. A stacked s pectrum of 11 
LBG s at z — 3.3 gave an upper limit of < 5% dBoutsia et al.l 
l201ll) . albeit dominated by a couple of bright LBGs, which may 
not be representative for the average high-redshif t galaxy popula- 
tion. lShaplev et al.ld2006b andlBoutsia et all j2oTlh use Starburst99 
templates i Leitherer et al J 1999b . with an intrinsic flux ratio of three 
between 1500 and 900 A. A larger flux ratio is the result of an older 
stellar population, which in turn depends on the assumed star for- 
mation history. If we assume the same intrinsic ratio (i.e. 3), the 
corresponding limit on the escape fraction derived for the M203 1 
galaxy would be a factor of two higher, i.e., / CS c,ll <22%. 

Even though lensing significantly boosts the signal from faint 
galaxies such that we can detect very small fluxes at rest-frame 
900 A for the M203 1 source, which has an intrinsic UV luminos- 
ity of 9L* determined from the SED fits in Paper I, we may have 
to look for intrinsically much fainter galaxies (<0.1L* ) that are 
expected to dominate the escaping UV flux at high redshifts. Fu- 
ture observations of gravitationally lensed faint galaxies offer an 
opportunity for determining / OS c,ll ( Vanzella et a lj2012l) . 
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Table 4. Physical properties of sources with emission lines 



Arc 


A31.1 


M0304 


M2031 


Section 


Magnification a 


26.6 ± 3.1 


42.0 ± 8.0 


15.8 ± 7.0 


Paper I 


2cm - 2 abs (km S _1 ) 


95 ± 131 


-10 ± 55 


80 ±90 


Paper I 


E(B - V) gas (mag) 


0.37 ±0.28 


0.23 ±0.01 


0.03 ± 0.33 


ED 


/esc,LL 


— 


— 


<11% 


EH 


/ csc,Lya 


85% 


0% 


20 - 50% 


□ 




21600 + 3000 


1 2900 + 900 


16100 + 300 


mi 


r) fO TTl /rm - ^ 
»tg \\J 11] ^LIIl ) 


144 ± 122 


Ol HI O 


QQn -L Ian 

OOU HI lOU 




n [C TTTl frm - ^ 
/ tg L>_ in J ^Liii J 










log i?23 


0.85 ± 0.05 


0.94 ± 0.001 


0.91 ± 0.02 


|3.4.2| 


12+log(0/H) (R 23 ) 


7.63 ±0.10 


8.18 ± 0.05 


7.74 ±0.03 


Paper I 


Z / Z @,R 2 3 


0.1 


0.3 


0.1 


Paper I 


12+log(0/H) (03JV2) 




8.04 ±0.01 




Paper I 


12+log(0/H) (AW02) 


7.85 ±0.15 


8.20 ±0.01 


7.56 ± 0.11 


Paper I 


12+log(0/H) (direct) 


7.69 ±0.13 


8.07 ±0.09 


7.76 ± 0.03 


|3.4.1| 


Z/Z®, direct 


0.07 


0.2 


0.1 




logtf 


-2.1 


-2.8 


-2.1 


|3.4.3| 


log(C/0) 


-1.03 ±0.08 




-0.80 ±0.09 


|3.4.3| 


log(N/0) 




-1.64 ±0.05 




13.4.41 



4 MODELLING THE Lya EMISSION LINES 

To infer the physical conditions governing the two systems for 
which Lya emission is detected, we undertake a series radiative 
transfer (RT) simulations, with the aim of fitting synthetic spec- 
tra to the observed ones. For this purpose, we appl y the three- 
dimen sional Lya scattering code MoCaLaTA i Laursen et ahl 
l2009ah . including the effects of dust dLaursen et alj|2009bl) . In the 
following, the basics of the simulations are outlined. 

The galaxies are modelled as spherical conglomerations of 
gas existing in two phases: warm, chiefly neutral clouds of rel- 
atively high neutral hydrogen density, dispersed in a hotter and 
more ionised intercloud medium (ICM). The ICM density, as well 
as the number density of clouds, decreases exponentially with dis- 
tance from the center. Additionally, the systems are surrounded by 
a spherical shell of outflowing gas. 

Densities, temperatures, and velocity fields are assigned to 
cells in a Cartesian grid. The grid has a base resolution of 128 3 
cells, with cells constituting the edge of clouds being refined re- 
cursively into eight cells, to make clouds more spherical. In addi- 
tion, Lya photons are also emitted from the centre and out with 
a luminosity which decreases exponentially with radius. A pho- 
ton is traced as it scatters its way out of the galaxy on individual 
hydrogen atoms, constantly changing direction and frequency, or, 
alternatively, is absorbed by dust. To obtain good statistics for the 
spectrum, a total of ~ 10 5 photons are emitted. 

At the redshifts of the two emitters (z — 1.8 and z — 3.5 
for A3 1.1 and M2031, respectively), the IGM is ionised to a high 
degree. Nevertheless, a fraction of the photons blueward of the Lya 
line may still be scattered out of the line of sight by diffuse neutral 
gas clouds in the circumgal actic medium. This e ffect is modelled 
in a statistical way following lLaursen et ai]( l201lh . 

In addition to matching the observed spectra, the preferred 
model should also be consistent with the star formation rates 
(SFRs) obtained in Paper I, as well as with the escape fractions 
inferred from comparing the integrated Lya and Ha fluxes. Under 
these constraints, and considering the large number of free param- 



eters in the model, rather than searching the full parameter space as 
has been done previously for a shell-only model of four parameters, 
including gas and dust column density, temperature, an d expansion 
velocity (Verha mme et alj |2008: Scha erer et alj|20l"lh . a series of 
simulations is conducted, changing the various input parameters 
until a satisfactory fit is obtained. This approach has been demon- 
strated to be able to matc h simultaneously a broad range of observ- 
ables for a single system dNoterdaeme et alj 2012). 

As an "initial" model, we make use of the few observational 
constraints we have at hand. Since we have measured the metallic- 
ities of the two systems (0.07 and 0.10 solar for A31.1 and M2031, 
respectively), this parameter is held fixed in the simulation to limit 
the already high number of free parameters. 

Assuming an SMC extinction law and a dust-to-metal ratio 
similar to the local universe (which is thought to be approximately 
valid even at high redshif t; see, e.g.. |Pei et alj [T999: Petti ni et"al] 
ll997l ; ISavagiioetalj2 003), the dust density is calculated. The Mo- 
CaLaTA code can model the dust as either LMC or SMC dust. We 
choose SMC dust, as the young stellar population of the SMC may 
be expected to match more closely those of the observed galaxies, 
but note that the difference between various extinction laws is of 
minor importance; over the narrow spectral region of the Lya line, 
th e dust cross section is almost flat, and indeed a comparison study 
in lLaursen et al. (2009b) of the two dust models reveals a change in 
escape fractions at the ~1% level only. Furthermore, the fact that 
the blue components of the spectra are significantly smaller than 
the red ones, albeit still visible, is a hint of modest outflow veloc- 
ities. Finally, for the temperatures we assume typical values of the 
three phases in the model: 10 4 , 2 x 10 4 , and 10 5 K for the clouds, 
the ICM, and the expanding shell, respectively. 



4.1 Best fit models 

Figure [2] shows the best fit models to the observed spectra, given 
the constraints. The synthetic spectra have been scaled by an SFR 
of 0.8 and 4.5 Mq yr" 1 , respectively, which is consistent with the 
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Figure 2. Observed (black) and fitted (red) spectra of the systems A3 1.1 (left) and M2031 (right). The yellow region shows the 68% confidence interval due 
to absorption in the IGM. 



Table 5. Galaxy model values 



Parameter 


A31.1 


M2031 


Galaxy radius 


1 kpc 


1 kpc 


Number of clouds 


80 


200 


Cloud radii 


20 pc 


20 pc 


Cloud temperature 


10 4 K 


10 4 K 


Cloud density" 


2.5 cm" 3 


100 cm -3 


Metallicity 


0.07 Zq 


0.1 Zq 


Cloud velocity dispersion 


50 km s _1 


10 kms" 1 


Central ICM density 


7.5xl0~ 2 cm -3 


3xl0~ 4 cm" 3 


ICM density" scale length 


0.4 kpc 


1 kpc 


ICM temperature 


5 x 10 4 K 


5 x 10 4 K 


Shell radius 


1 kpc 


1 kpc 


Shell expansion velocity 


60 km s _1 


110 kms" 1 


Shell column density" 


3 X 10 18 cm -2 


10 18 cm" 2 


Shell temperature 


5 x 10 s K 


5 X 10 s K 


Lya emission scale length 


0.1 kpc 


0.25 kpc 


Intrinsic Lya EW 


40 A 


20 A 



"All densities refer to neutral hydrogen only. 



values found in Sect. [3] Moreover, the escape fractions were found 
to be 0.87 and 0.56, respectively. While the former matches per- 
fectly the value found for A3 1 . 1 , the latter is somewhat higher than 
the ~0.20 found for M2031, although due to the rather large uncer- 
tainty on the fluxes, it is not inconsistent. Table [5] summarises the 
final parameter values used for the models. 

Note that the models are perfectly scalable: radii of 1 kpc and 
2 kpc have been assumed for A3 1 . 1 and M203 1 , respectively; how- 
ever, the exact same spectra and escape fractions would emerge 
from models with twice the radius, but with all densities a factor of 
4 smaller, and scale lengths and cloud radii doubled. 

As mentioned in Sec. [2] the asymmetric profile is characteris- 
tic of an expanding ISM. However, the fact that we do see the blue 
peaks in both spectra tells us that the typical outflow velocities must 
be quite low. For a homogeneous medium, the blue peak vanishes 
already at V ou t ~ 100 kms" for a broad range of column densi- 
ties (e.g. [Verhamme e t alj|20060 . For an inhomogeneous medium, 
since the photons blueward of the line may find low-density paths 



out of the galaxy, higher velocities may be allowed with the blue 
peak still visible. 

The outflows of the two models that match the observed spec- 
tra best are 60 and 110km s ~ 1 forA31.1 and M203 1 , respectively. 
This is consistent with the values inferred from comparing the red- 
shifts of the absorption and emission lines (see Tab.[4] Although the 
rather large uncertainties in the latter values make them consistent 
with zero as well, we consider it likely that the two emitters in fact 
do exhibit outflow. In principle, foreground diffuse H I clouds could 
also cause a decreased blue peak. At the given redshifts the IGM 
is likely to be highly ionised. However, its clumpy nature could in 
principle cause absorption nevertheless. In Fig. [2] the yellow region 
marks the 68% confidence level of the Lya emission line, after be- 
ing transmitted through the IGM. Indeed, the expected impact of 
the IGM is rather low. 

Most observed Lya emitters exhibit an asymmetric line pro- 
file with the blue peak missing. However, with sufficient spectral 
resolution it seems that the fraction of Lya e mitters where the 
blue peak is still vi sible is quite large (20-50% iKulas et aljkod ; 
lYamada et alj2012b . 

It is important to keep in mind that we do not claim that the 
model galaxies must be a realistic representation of the real galax- 
ies. Many processes shape the Lya line, and these are to some ex- 
tend degenerate. Nevertheless, it is striking that it is possible to 
match the lines under the restrictions set by the other observables 
such as metallicity and SFR. 

The higher SFR and metallicity of M2031 relative to A3 1.1 
hints at the former being more evolved than the latter. This inter- 
pretation is backed up by the results of the simulations; the higher 
number of clouds and the larger density contrast between clouds 
and ICM implies a more fragmented ISM. The shallower density 
gradient could then be caused by the enhanced star formation "in- 
flating" the galaxy, while the higher wind velocity is a natural con- 
sequence of a more vigorous feedback. 

5 SUMMARY 

We have measured the oxygen abundance using direct T e methods 
for three galaxies at z = 1.8, 2.0 and 3.5, respectively. The direct 
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oxygen abundance measurements are in agreement with those de- 
rived f rom the conventional R23 ratio using the diagnostics calibra- 
tion in IPilvugin & Thuarj d2005h . This proves that the strong-line 
diagnostics are in fact useful to determine oxygen abundances in 
high-redshift galaxies, as the temperature-sensitive lines for field 
galaxies are almost always below the detection limit even on 8-m 
telescopes. 

In addition, we have measured the C/O and N/O ratios for the 
three galaxies, and argue that their ratios are consistent with the 
general trend of sub-solar metal ratios for low-metallicity galaxies, 
and in particular for young galaxies where contributions from either 
nitrogen or carbon from older stars have not yet occured. 

We have modeled the emission line profiles in two strong Lya 
emitting galaxies using the constraints (SFR, reddening, escape 
fractions and metallicity) obtained from the analysis of their spec- 
tra. The relatively small inferred expansion velocities of the neutral 
gas shells correspond to the small velocity offsets seen in absorp- 
tion and rest-frame emission lines. The observed absorption and 
emission line redshifts show that the velocity offsets are consistent 
with zero at the la confidence level, but are also consistent with the 
small shell expansion velocities (60 and 110 km s" 1 , respectively) 
found in the radiative transfer simulations. Compared to more mas- 
sive LBGs with higher SFRs we observe a smaller velocity offset 
between the absorption and Lya emission wavelengths. To develop 
our understanding how galaxy outflows are related to the SFRs and 
escape of Lya photons, it would be interesting to expand the sample 
and verify if the trend of dec reasing velocity offs et with increasing 
Lya EW observed for LBGs JShaplev et al.l2003h is present also in 
Lya emitters with much higher EWs. 

Recent studies of Lya emitters have shown that the absorp- 
tion lines are more difficult to detect at progressively larger Lya 
equivalent widths, even when stacking the data from many emitters 
dBerrv et alj|2012l) . Gravitational lensing is the only means to ob- 
tain med ium to high-resolution spectroscopic data of high-EW Lya 
emitters jDiehl et alj2009MOuider et alj2009l ; rBavliss et alj201Ch . 
The A3 1 . 1 source, with an equivalent width of 40 A, does not ap- 
pear to have strong absorption lines, but the faint continuum and 
low signal-to-noise ratio of the spectrum prevents a more detailed 
analysis besides the strongest Sill lines at 1260, 1264 and 1526 A. 
Since one hour integration time per target was used in this study, 
we conclude that it would be possible to investigate the detailed 
physics of a strong Lya emitter with a few more hours integration 
time on a large telescope. 

Future observati ons of strong lensin g clusters, such as the 
HST/CLASH survey jPostman et alj|2012h . may reveal more po- 
tentially interesting lensed high-redshift sources with steep UV 
slopes and high surface-brightness regions. These could be the best 
targets to locate and investigate the interstellar absorption lines in 
strong Lya emitting galaxies. These galaxies are also the prime 
candidates to find other strong UV lines such that abundances of 
several species could be determined by direct methods. With such 
observations, we could constrain better chemical evolution models 
of high-redshift galaxies. 
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